Superconducting pairing and density-wave instabilities in quasi-one-dimensional 
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Using a renormalization group approach, we determine the phase diagram of an extended quasi- 
one-dimensional electron gas model that includes interchain hopping, nesting deviations and both 
intrachain and interchain repulsive interactions, d-wave superconductivity, which dominates over 
the spin-density-wave (SDW) phase at large nesting deviations, becomes unstable to the benefit 
of a triplet /-wave phase for a weak repulsive interchain backscattering term gi > 0, despite 
the persistence of dominant SDW correlations in the normal state. Antiferromagnetism becomes 
unstable against the formation of a charge-density-wave state when gi exceeds some critical value. 
While these features persist when both Umklapp processes and interchain forward scattering (gjt ) 
are taken into account, the effect of g 2 alone is found to frustrate nearest-neighbor interchain d- 
and /-wave pairing and instead favor next-nearest-neighbor interchain singlet or triplet pairing. We 
argue that the close proximity of SDW and charge-density-wave phases, singlet d-wave and triplet 
/-wave superconducting phases in the theoretical phase diagram provides a possible explanation for 
recent puzzling experimental findings in the Bechgaard salts, including the coexistence of SDW and 
charge-density- wave phases and the possibility of a triplet pairing in the superconducting phase. 

PACS numbers: 71.10.Li,74.20.Mn,74.70.Kn 



I. INTRODUCTION 

The theory of low-dimensional metals has exerted a 
strong influence on our understanding of ordered phases 
in quasi-one-dimensional (quasi-lD) organic conductors. 
The description of the extensively studied Bechgaard 
salts series ((TMTSF) 2 X) and their sulfur analogs, the 
Fabre ((TMTTF) 2 X) salts, has served over more than 
two decades to illustrate this viewiiS^ For these mate- 
rials, a direct correspondence can be traced between the 
various modulated spin and charge ordered states of their 
phase diagram and the possible states of the quasi- ID 
electron gas model when the couplings are repulsive and 
the filling of the band is commensurate with the underly- 
ing latticeA^ In its standard form, this generic model is 
defined by the phenomenological expression of the direct 
interaction between electrons in terms of weak intrachain 
backward (<?i) and forward (52) electron-electron scatter- 
ing processes, to which Umklapp ((73) scattering ampli- 
tudes are added at commensurate band filling. 5 ' 6 ' 7,8 The 
quasi- ID character of the model is defined by adding an 
interchain single electron hopping integral tj_ , which is at 
least an order of magnitude smaller than its longitudinal 
counterpart. 

In virtually all compounds of the above series, su- 
perconductivity is observed next to a spin-density-wave 
state for some critical value on the pressure scale(2*i2*ii*i2i 
whereas antiferromagnetic spin correlations are found to 
dominate the metallic state precursor to superconductiv- 
ity over a wide range of temperatures pii^ii On the theo- 
retical side, however, for spin-independent repulsive cou- 
plings and for a Fermi surface with good nesting prop- 



erties, the coexistence of spin-density-wave (SDW) and 
superconducting (SC) correlations is essentially excluded 
from the model phase diagram^is It is only when devi- 
ations from perfect electron-hole symmetry (nesting) are 
introduced and the long-range component of the SDW or- 
der is suppressed that superconductivity can be actually 
realized in place of magnetismii 7 ^^ In the anisotropic 
metallic phase where this suppression takes place, in- 
terchain Cooper pairing is enhanced and superconduc- 
tivity emerges from the coupling between the weakened 
electron-hole and the still singular electron-electron scat- 
tering channels. Recent calculationsi2i22*2I using the 
renormalization group (RG) method did confirm the exis- 
tence of such an electronic pairing mechanism beyond the 
level of single-channel RPA-like approachesi 22 i 23 i 24 i 25 i 26 
A smooth crossover from the SDW state to supercon- 
ductivity has then been found whenever the amplitude 
of nesting deviations reaches some threshold - a result in 
accordance with the sequence of transitions observed as 
a function of pressure AiLSJ. 

The RG approach shows nevertheless that whenever 
antiferromagnetism stands out as the dominant corre- 
lation in the normal state, the most stable interchain 
pairing is invariably a spin singlet state corresponding 
to a 'd-wave' symmetry gap with nodes on the Fermi 
surface. Although some experimental findings in the 
Bechgaard salts and their sulfur analogs do agree with 
this type of pairing, 2 L2£. other series of observations have 
rather been interpreted in support of a triplet order 
parameter ^S^SL^ thus challenging the singlet scenario 
for superconductivity. Though there is too little so far 
to favor one scenario over another, these observations 



bring us to the question of whether triplet superconduc- 
tivity can be possible or not when short-range antiferro- 
magnetic correlations are dominant in the metallic state. 
As is well known for the electron gas model, there is a 
region of the phase diagram where triplet 'p-wave' su- 
perconductivity does exist as the most stable state. It 
has been suggested, on a phenomenological basis, that 
such a state is realized in the Bechgaard salts, 33 Within 
a microscopic approach, the region where p-wave super- 
conductivity is stable is defined by irrelevant Umklapp 
scattering and by a backward scattering coupling that is 
much larger than forward scattering. However, this de- 
scription pattern for superconductivity can be considered 
unsatisfactory given the unrealistic constraint it puts on 
the coupling constants, and for the suppression of both 
the antiferromagnetism at short distance and the Mott 
pseudo gap in the charge sector^ These flaws can hardly 
be reconciled with the related phenomenology of these 
molecular compounds observed around the critical pres- 
sure for superconductivity, when either the temperature 
or the magnetic field is variedi 1 ^ 4 * 3 ^* 3 ^ 

Another way to look at this problem is to consider 
more closely the effect of charge fluctuations on super- 
conductivity. Staggered charge fluctuations are known 
to favor triplet pairing at odd but large angular mo- 
mentum, a mechanism whose roots go back to the early 
work of Kohn and Luttinger about Cooper pairing in 
the presence of charge - Friedel - oscillations in isotropic 
Fermi systems^ For the quasi-lD electron gas model and 
its version for lattice electrons, recent calculations have 
shown that triplet /-wave pairing is indeed enhanced 
when intrachain couplings are chosen to boost charge- 
density-wave (CDW) fluctuations close to the level found 
in the SDW channeliSi* 3 ^ 3 ^*^ However, for realistic re- 
pulsive couplings, d-wave pairing still remains tied to the 
highest critical temperature and hence to the most stable 
state for superconductivity. 

All this goes to establish the robustness of d-wave pair- 
ing for the model with repulsive intrachain couplings and 
nesting deviations. It turns out, however, that the model 
is incomplete when charge fluctuations are found along 
the chains since then interchain Coulomb interaction is 
also present in practice. The inclusion of direct interchain 
electron-electron scattering processes, which will be de- 
noted by 3^1.2,3 m the following, defines the quasi-lD 
electron gas model in its extended form4i At large mo- 
mentum transfer, the interchain interaction is well known 
to favor a CDW ordered statei 41 i 42 i 43 i 44 ' 45 This mecha- 
nism is mostly responsible for CDW long-range order ob- 
served in several organic and inorganic low-dimensional 
solids. 4 £ji£ The physical relevance of interchain interac- 
tions in the Bechgaard salts, besides the intrachain gi and 
t_i_, is supported by x-ray studies, which revealed that the 
SDW phase of these compounds is actually accompanied 
by CDW order4242iii2i On the theoretical grounds, very 
little is known about the impact of adding direct inter- 
chain interactions on the structure of the phase diagram, 
especially in the repulsive sector when both a finite t± 



and nesting deviations are present. 

In this work we wish to determine the possible density- 
wave and superconducting states of the extended quasi- 
lD electron-gas model. In order to tackle this problem 
we shall apply the renormalization group method, which 
at the one-loop level has proved to be suited to reach 
a controlled description of interfering density-wave and 
superconducting channels of correlations. Among the re- 
sults reported below, we have the unexpected finding that 
a small repulsive interchain backscattering term g^ > 
is sufficient to make d-wave superconductivity unstable 
to the benefit of a triplet /-wave phase. This occurs de- 
spite dominant SDW correlations in the metallic state 
and stable itinerant antiferromagnetism at lower nest- 
ing deviations. Under the latter conditions, SDW order 
becomes in turn unstable to the formation of a CDW 
state when the amplitude of g^ exceeds some critical 
value. While these features persist when commensur abil- 
ity effects are taken into account and small - half-filling 
- Umklapp scattering is included, the effect of interchain 
forward scattering (g^ > 0) is found to frustrate nearest- 
neighbor interchain d- and /-wave pairing and to favor 
instead superconductivity with next-nearest-neighbor in- 
terchain pairing. Part of these results have been reported 
in Refs. lslHl 

In Section [HJ we introduce the model and the RG 
scheme employed for the four point vertices and the re- 
sponse functions. By way of illustration, the RG results 
at the one-loop level are given for purely intrachain inter- 
actions. In Section lTTTl we present the results for non zero 
interchain backward and forward interactions by which 
the different possibilities of ordered states in the phase 
diagram are obtained in the incommensurate case. The 
influence of Umklapp processes in the half-filled case is 
examined in Section IIVI A discussion of the results is 
given in section [3 where a possible connection between 
theory and experiments is made. 



II. THE EXTENDED ELECTRON GAS MODEL 
A. Model 

We consider a lattice of N± coupled metallic 
chains described by the partition function Z = 
JJ Dip* Dip e~( So+Sl \ which is expressed as a functional 
integration over anticommuting (Grassmann) fields ip. 
Here So and Si are the noninteracting and interacting 
parts of the action, respectively. The former part is given 
by 

So = - ^2^K( ik - £.k)lpK , (1) 
K 

where K — (a, k) , k — (fco;k), ko is the fermionic Mat- 
subara frequency, k = (kn,k±), the wave vector, and 
(7 = ±, the spin of the fermion field. The kinetic energy 
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is given by 

6c = e-k- fi = v F (\k\\\ - k F ) + ex(kx) , (2) 
e±(k±) = — 2t± cos k± — 2t' ± cos 2k± , (3) 

where /i is the chemical potential taken as temperature 
independent. Here we have linearized the spectrum in 
the chain direction, using vp = 2tii sinfci? as the longi- 
tudinal Fermi velocity, kp being the parallel Fermi wave 
vector if t± = 0. Throughout this work both the chain 
lattice constant a and the interchain distance d± are put 
equal to unity. The interchain single electron hopping ij_ 
is considered small with respect to the longitudinal band- 
width 2Aq. In the following, we will take 2Ao = 30tj_, 
which is a typical figure for the anisotropy ratio in quasi- 
1D conductors like the Bechgaard and Fabre saltsS^ 
The next nearest-neighbor hopping t' ± <C ij_ in the trans- 
verse direction gives the amplitude of nesting deviations. 
In Eq. J3J), we have neglected the possibility of hopping 
in the third direction, which does not have any sizeable 
effect on our calculations. Its existence is of course cru- 
cial for the stabilization of true long-range order at finite 
temperature. 

In the framework of the quasi-lD electron gas 
model)£*i& the electron-electron interaction is parameter- 
ized by means of the g-ology approach. One first dis- 
tinguishes between right and left moving fermions, de- 
pending on their velocity along the chains, so that the 
Grassmann variables become 

Using these definitions, the interaction part of the action 
takes the form 

Si = J-f ^2 y^M+fcj.fca+fci mod G (5) 
k' 1 k' 2 k2ki <j<j' 

x^g^k'^k^Rl^Ll^Rk^-Lk^ 
+92 (k[ k' 2 k 2 h) R* k{ a L* k , 2<T ,L k2(T ,R kl(T 

+\{93{k[k' 2 k2ki)Rt fia Rl^,L k2a ,L kia + c.c.)| , 

where N is the number of lattice sites and G = (0; 4fcp, 0) 
the reciprocal lattice vector along the chains, which is in- 
volved in the Umklapp scattering at half-filling. Here we 
have omitted the so-called 34 contribution for the inter- 
action of electrons of the same branch since this coupling 
does not contribute to the singular channels of correla- 
tion in the renormalization group flow at the one-loop 
level. If we restrict ourselves to intrachain and nearest- 
neighbor chain interactions, the amplitudes of the bare 
interactions are given by 

gj(k' ±1 k' ±2 k ±2 k±i) = gj + 2gf cos(k' xl - feu) , (6) 

where the only - transverse - momentum dependence 
comes from the interchain interaction. In this work we 



analyze the properties of the model for repulsive gj , gj~ 
and finite Umklapp scattering, which is the physically 
relevant sector for real materials like the Bechgaard salts 
and their sulfur analogs. 



B. Renormalization group equations for the 
interactions 

The renormalization group can be applied to the prob- 
lem of quasi-lD systems of interacting electrons in dif- 
ferent ways>i2iSi££i££ In the following, we shall adopt 
the so-called one-particle irreducible (1PI) momentum 
shell RG scheme, as developed by Honerkamp et al*& (cf. 
also Binz et al£2>). In this scheme, the 1PI vertex func- 
tions T\(Kj) of a physical system with infrared cutoff 
A are calculated. Only degrees of freedom with energies 
equal to or greater than the cutoff are integrated out, 
which allows one to establish a direct link with the low 
energy effective action Srjf with an ultraviolet cutoff, 
as obtained in a Kadanoff- Wilson renormalization group 
scheme. According to Morris, 60 if ekJ < A, the effec- 
tive interactions g e ^ (Kj) that appear in S e ^ are simply 
g e A ff (K j )=r A (K j ). 

The free particle propagator is suppressed for energies 
below an infrared cutoff A, 

Ck = {Mtc)s = ©(I6c| - A) x — 1— , (7) 

iko - Ck 

where O is the step function. Setting all vertices involv- 
ing more than two particles to zero, the renormalization 
group equations to the one-loop level are obtained. We 
also neglect self-energy corrections. 

It is practically impossible to take into account the 
complete functional dependence of the vertices, which 
will in general vary with the wave vectors k as well as 
with the frequencies fco of the incoming and outgoing par- 
ticles. By means of scaling arguments, one can show that 
the frequency dependence is irrelevant in the renormal- 
ization group sense, cf. for example Ref. l6ll The same 
is true for the dependence on the distance of the wave 
vector k from the Fermi surface. Irrelevant variables are 
not necessarily negligible, in the sense that they may in 
principle have an influence on the RG flow before they 
vanish. However, the k s dependence, with k s the projec- 
tion of the wave vector on the Fermi surface, remains by 
far the most important one in the low energy limit. We 
will hence only consider the dependence of the functions 
r on the positions of the wave vectors along the Fermi 
surface, which we will parameterize by k±. Due to mo- 
mentum conservation, it is in general impossible to have 
all four arguments of a given vertex on the Fermi surface. 
Usually, the flow of the vertices with three arguments on 
the Fermi surface is calculated. In the case of imperfect 
nesting, such a procedure will underestimate the SDW 
correlations. In our case, the strong anisotropy of the 
dispersion relation allows us to make a slightly different 
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choice. We fix the values of k\\ and kg such that the most 
divergent contributions to the renormalization group flow 
are taken into account. If we define the Cooper, Landau, 
and Pcicrls momentum variables 



(8) 
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we can write the one-loop RG equations in the form 



fcloo 



Peier\s(T(qc,qc,qv)) 



In this equation, qc, qc and q-p are functions of the in- 
ternal loop momentum- frequency k\ oop . The most di- 
vergent contributions to the RG flow are obtained if 
qc, o — q^-po — for the frequencies, as well as 

qc\\ = , q v \\ = 2k F for r\ 2 , (9) 
q!p\ = ±2k F for T 3 , 



for the longitudinal components. Note that, for either 
choice of the k» dependence, one does not obtain a closed 
set of equations. The renormalization group equations 
for vertices for which relations © are valid contain ver- 
tex functions with different values of qc, q-p- We will 
therefore need to replace these functions by those calcu- 
lated by imposing conditions JjJJl. 

For the dimensionless vertices T = T/ttvf, we thus 
obtain the equations^ 



Ti{k' xl k' X2 k X2 k xl ) 



k± 

X {f l(V xl V X3 k±V ± )f 2(k X k' x k X1 k X2 ) 
+f 2 (k xl k X2 k' x k X )f 1 (k X k x k X2 k X1 )} 

-E" 



k' ± = -k ± +q c± 
1C± =k ±l+ k ±2 



B v (k x ,q vx ) 



x{(f 2 (A;i 1 fcxfcxifci)-fi(^ 1 fcxA;lA;xi))f 1 (A;lfci 2 fe ±2 A; ± ) 

+f 1 (k' xl k x k' x k xl )(f 2 (k , x k' X2 k x k X2 ) - f 1 (fc ± fc ±2 A;x2fc±)) 
+ (T 3 (k' xl k x k xl k' x ) - f 3 (k' xl k x k' x k xl ))T 3 (k' x k' X2 k X2 k x ) 

+f 3 (k' xl k x k' x k Xl ) (f 3 (k' x k' X2 k x k X2 ) - f 3 (k' x k' ±2 k X2 k x ))} 



(10) 



ir'n 



T 2 {k' xl k' X2 k X2 k xl ) 



B c(k±,qc±) 

k± L 

*{T 1 {k' xl k' X2 k x k' x )f l {k x k' x k xl k X2 ) 

+r 2 (k' xl k X2 k x k x )f 2 (k x k x k X2 k xl )} 

+ TrYl B v {k x ,q vx ) 

x [r 2 (k' xl k x k X2 k' x )f 2 (k' x k' X2 k x k X i) 

+f 3 {k' xl k X k X2 k' x )f 3 (k' x k' X2 k X k X1 )} 



k ± =-k ± +q c± 
1C±= k -Ll+ k -L2 



k ± =k ± +q-p ± 
1-P±= k '±l~ k ±2 



(11) 
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T 3 (k'_ L1 k'_ L2 k± 2 k±i) 



-E 



{^2(k' ±1 k±kxik' ± ) - f i(fcxi^xfcxfcxi))f 3(k' ± k' ±2 k±2kx) 
+Ti(k'j_ 1 k^k' J _k±i)(r 3 (k' ± k' ±2 kj_kj_2) - T 3 (k' ± k' ±2 kx2kx)) 
+ {f 3 (k' ±1 k ± k ±1 k' ± ) - t 3 (kf ±1 k x k^k ±1 ))r 1 ^ ± k f X2 kr i kx) 

+f 3 (^ 1 fe ± fc^fc ±1 )(f 2 (fcl^ 2 fc ± fc ±2 )-f 1 (^fel 2 fc ±2 fc±))} 
- 1 k± L 

x {f 2{k' ±1 k ± k X2 k' ± )f 3 (k' ± k' ±2 k x k ±1 ) 
+f 3 (k' ±1 k ± k± 2 k' ± )f 2 (k' ± k' ±2 k ± k ±1 )\ f , 

q v± =k ±1 -kj_ 2 



(12) 



k' ± =k 



where the dot denotes derivation with respect to 
— ln(A/A ). The particle-particle (C) and particle-hole 
(P) loops, after summation over Matsubara frequencies 
and longitudinal wave vectors, are given by 



B c/V (k±,q±) = E 0(| A + ^ 



c/v\ 

u=±l 

If A + vA cjv 

x - tanh ■ — 

2 V 2T 

A 



A) 



(13) 



u A \ 

tanh — 
2T J 



Ac(kj_,q±_) = 
A v (k_L,q±) = 



2A + vA 

-e x (k±) 
-e±(A±) 



ex(~k± + q± 
e±(k± + q±) 



(14) 



For continuity reasons, we always take 8(0) := -i. The 
starting values are given by the bare interactions, cf. sec- 
tion \nr\ 

It is useful to consider the following limiting cases for 
the RG equations. As long as A t±, we may take 
their ID limit (tj_ = 0). In this regime, the Cooper 
and Peierls renormalization channels are entirely coupled, 
in the sense that all vertices are strongly renormalizcd 
in both channels (except for Umklapp processes, which 
only appear in the Peierls channel). When A « ii, 
the coupling between the Cooper and Peierls channels 
is weak in the sense that, depending on the arguments 
k±, most vertices are strongly renormalized in only one 
(or no) channel at a time. Nevertheless, the remaining 
interplay between the channels is at the origin of spin- 
and charge-fluctuation-induced superconductivity in the 
weak coupling regime. For purely repulsive interactions, 
the Peierls channel is the most important one as long as 
deviations from perfect nesting may be neglected, i. e. 
for A t' ± . In the presence of attractive effective inter- 
actions, the Cooper channel plays an important role, and 



will be dominant when A < The RG equations are 
written in these limits in appendix 1X1 

A similar dimensional crossover may be observed with 
respect to the temperature instead of A. As long as 
T tx, the functions we calculate do not vary very 
strongly with the transverse wave vector, contrary to the 
low temperature case. An example is given in Fig.^a). 



C. Response functions 

In order to evaluate susceptibilities, we add a term of 
the form 

^ = EX>aSc(?)OaSc(5) +CC. (15) 
a q 

+ EE h ^)Oi M nUq)+c.c. 

a,M q 

to the action. The first term describes pairing and the 
second term density- wave correlations. The external 
fields h,(*\q) are taken to be infinitesimal. They cou- 
ple to pairs of fermionic variables which we will define in 
the following. Let us first introduce the particle-particle 
operators 



(16) 



for singlet (a = s) and triplet (a = t X} y :Z ) pairs. The 
spin dependence is given by the coefficients 



( s ) x 



(17) 



-(**) _ 



l *> - -aS i r {ty) - -iS i t kz *> - 5 



(*») _ 



(<7 = ±1). The pair operator appearing in Sh is defined 
as 

O a sc{q) = \jjj E z a (q - k,k)o a (k,q) . (18) 



{fc;fc,|>0} 
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name 


spin pairing 


A P (Ax) 


s 


singlet 


1 


Px 

Py 


triplet 


r 

sin k± 


dxy 


singlet 


cos k±_ 
r sin k± 


f 


triplet 


r cos k± 
sin zk± 


9 


singlet 


cos 2fcx 
r sin 2k± 


h 


triplet 


r cos 2fcx 


i 


singlet 


cos 3k± 



TABLE I: Superconducting order parameters in a quasi-lD 
geometry. The names are assigned according to the number 
of sign changes along the Fermi surface. 

The function z a describes the orbital symmetry of the 
superconducting order parameter. We classify these or- 
der parameters by their behavior on the Fermi surface. 
They are parameterized by k± and r = sgn k\\ . A list of 



with M = + for site and M = — for bond density 
waves, and Za (k,q) = 1. (For the case M = — 1, 
we have neglected a constant imaginary factor.) Note 
that (0i M) )* = MO ( a M) , and that (0 ( a M) ( a M) ) and 
(Oa^*0^*'*) do not vanish, but contribute to the as- 
sociated susceptibilities. Away from half filling, we only 
consider site density waves. In this case, only the first 
term in Eq. (1221) contributes. 

In the presence of Sh, the one-particle vertex con- 
tains a nondiagonal part with a contribution linear in 

( h aDw(l)T ° f the form 

^4^- 9 Ci M) (fc',fc) (23) 



the superconducting order parameters we shall examine 
is given in Table [I] With these, we have 

, v _ J 1 for s and p x , 

Similarly, we introduce particle-hole operators 

o a (k,q)=J2^lLt- q ,,'Rk tir (20) 

a' a 

for charge (a = C) and spin (a = S x ,y,z) excitations. 
Here, is the 2x2 identity matrix, and cr^" ^ are 
the Pauli matrices. At half filling, one has to distinguish 
between bond and site density waves. In direct space, 
these are given by 

O = /3E^^k,A' for site DW, 

X ||E^(^k^^ + c. C .) for bond DW, 

(21) 

where d = (0; d|| , 0), and d\\ is the lattice periodicity along 
the chains. Since we have kpd\\ = ^, the associated 
Fourier transforms are, for q» rss 21cf, 



(22) 



I 

(if k\i > 0). Also due to Sh, vertices with two outgoing, 
but no incoming, particles are now non-zero. Their linear 
part in (h a sc(q))* takes the form 



rj#** jg _ fc Ca(A',*) (24) 



(if fcii > 0). The vertex parts C Q determine the renor- 
malization of the susceptibilities Xa as follows. In the 
beginning of the flow, we have £ a (k',k) = z a (k' ,k) and 
Xa = 0. The RG equations for density waves (a = C, S) 
are: 



J 



oiVwM^l]^ E [z i ™Xk- q ,k)o a (k, q )+M(zi M \k + q -G,k))*o* a (k,G-q)] 

{fc;fc||>0} 



6 



&(p±,P±+q±) = ^VB P (^9i)( a M (fci,fci+«i) (25) 
x [T a (kx + qx,Px,k±,p± + q±) 

,p±,k± +q_L,Px + q±)] , 

For Cooper pairs (a = s,t), we have: 

( a (-p± + q±,P±) = -Ty-^2B c (k±,q±)( a (-k± +q±,k±) (27) 
X k± 

xf a (k±,-k± + q±, -p± + q±,p±) , 
2 



Xa(q±) = -^X)Bc(fc±,^)|Ca(-fc±+<^,M 2 . (28) 

k± 

X is defined as nv F ^x- We have introduced the linear combinations of vertex functions 

^c(k' ± i,k' ±2 ,k ±2 ,k ±1 ) = -2T 1 {k , 1 _ 1 ,k' ±2 ,k^ l ,k 1 _2)+Y 2 {k' ±1 ,k' ±2 ,k^2,k^ l ) , (29) 
T s (k' ±1 ,k'j_ 2 ,k ±2 ,k±i) = T 2 (k' ±1 ,k' ±2 ,k ±2 ,k±i) , 

r 3 C) ( fc ±l^±2> fc ±2, kj_l) = -2r 3 (k' xl ,k' X2 ,k±l,k_L2) +r 3 (k' ±1 ,k' X 2,k±2,k±i) , 
r 3 ( fc ll7^2; fc i2,fc±l) = r 3 (fc_ L1 ,fe_ L2 ,A;_L2,fc±l) , 

r s (k' ± i,k'j_ 2 ,k±2,k±i) = -T 1 (k' ±1 ,k' ± 2,kj_i,k±2) - T 2 (k' A _ 1 ,k' 1 _ 2 ,k_ L2 ,k±i) , 
r t (fc^ L1 ,A:^ L2 ,fc ± 2,fc±i) = Ti(k , ±1 ,k' ±2 ,k±i,k±2) — T2{k' 1 _i,k'±2,kAj2,k±i). 

Interactions g a (a — C,S X! y iZ ,s,t XtVvZ ) and g% ,g 3 are defined analogously. Using these functions as well as the 
particle-particle and particle-hole pair operators defined above, and neglecting Umklapp processes, we can rewrite the 
interaction part of the action in the following way: 

Si = -- ^2 J^^29a(k',k- q,k' - q, k)o* a (k' , q)o a (k, q) (30) 

ot=C,S x y jZ k' kq 



= ~9 X] at ^29a(k',q~ k',q- k,k)o* a (k' ,q)o a (k,q) , 



2 V N 



g St t (resp. gc,s) thus describes the interaction between 
particles forming a singlet or triplet pair (resp. particle- 
hole pair). The relation between these couplings is as 
follows: 

9s = ^(3gs-gc) , (31) 
9t = \{gs + gc) ■ 

We will make use of these relations in the following, 
when we discuss density fluctuation induced Cooper pair- 
ing. Note that Umklapp processes do not couple to the 
Cooper channel, and therefore affect Cooper pair forma- 
tion only indirectly via their effect on gc, gs- 

The RG equations are solved numerically, using a 
fourth order Runge-Kutta algorithm^ 3 , with fixed step 



sizes (taking large steps in the beginning of the flow and 
short steps close to the divergence). The graphs show- 
ing transition temperatures are obtained with the help 
of an adaptive stepsize algorithmiSiSi As far as the k± 
dependence of the vertex functions is concerned, we dis- 
cretize the Fermi surface using 32 patches for each sheet. 
Taking advantage of all symmetries of the problem, we 
thus have to calculate 9010 different function values for 
each Tj . The use of twice the number of patches does not 
significantly modify our results. The integral of the func- 
tions Bq,v must be calculated more precisely; we use a 
fourth order Runge-Kutta algorithm with adaptive step 
size. 
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FIG. 1: RG phase diagram for the quasi-lD electron gas 
model at incommensurate filling (gz — 0) with intrachain in- 
teractions (gf = 0) and nesting deviations, parameterized by 
t'j_. The circles indicate the transition temperature for SDW 
and the triangles that for SCd. All figures in this article 
are obtained using the bare intrachain interactions g\ = 0.32, 
<?2 = 0.64 and the anisotropy ratio Ao/tx = 15. 



D. Results for the case of intrachain interactions 

To illustrate our method, we shall first consider the 
by now well known case of a quasi-lD electron gas model 
with purely intrachain interactions and no Umklapp scat- 
tering. This case has been studied in detail by Duprat 
and BourbonnaisA 9 - using a Kadanoff- Wilson RG scheme 
and an approximation where only two independent mo- 
mentum variables for the vertices were taken into ac- 
count. We here use the 1PI scheme and retain the full 
three variable dependence for the vertices. Our results 
confirm those of reference fl9l 

The RG calculations that follow are performed for the 
values gi = 0.32 and g2 = 0.64, which are representa- 
tive of the couplings likely to be found in practice in 
low dimensional conductors. Thus for not too large nest- 
ing deviations, the renormalization group flow scales to 
strong coupling, which leads to a singular behavior in 
the susceptibility of a particular channel of correlations. 
This signals an instability of the normal state towards 
an ordered phase. The phase diagram obtained as a 
function of t' ± is shown in Fig. ^ For good nesting, 
there is an SDW phase with a modulation wave vec- 
tor Qo = (2kp,ir) that corresponds to the best nesting 
vector of the spectrum The transition temperature 
obtained for perfect nesting (t' ± = 0) is T c ° w 0.055ix- 
If we take t± « 200 K, we have T c ° w 11 K, which falls 
in the range of the experimental T c for systems like 
the Bechgaard salts at low pressure. By increasing t' ± , 
the transition temperature decreases until the threshold 
value t' Xc w 0.045£j_ ~ 9 K is reached, where the SDW 
is suppressed and replaced by d-wave superconductivity 
(SCd). The maximum temperature for the SCd state 
is found to be T c °(SCd) w 0.002ix w 0.4 K. These esti- 
mates for t'j_ and T c are comparable to the experimental 
results, 9 i 19 i 65 i 66 
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FIG. 2: Temperature dependence of the susceptibilities in the 
normal phase (gz = 0, gf = 0): 

(a) above the SDW phase (t'± — 0.04tj_). The continuous 
(dashed) line corresponds to the SDW response at q± — n (0); 
the dashed-dotted line to SCd correlations. It is interesting 
to observe that the temperature scale T ~ t± below which 
the q± — curve separates and levels off corresponds to the 
so-called single particle dimensionality crossover i— Our RG 
scheme thus captures this effect correctly. 

(b) above the SC phase (t'± = 0.048t±). The continuous 
line corresponds to the SDW response with transverse mod- 
ulation 7r, the dashed-dotted line to SCd a; 2_ !/ 2 (cosfcx), the 
dotted line to SCd xy (r sinfcx) and the dashed line to SCg 
(r sin2fcx) correlations. 



The order parameter of the low temperature phase can 
be identified in two different ways, which give equivalent 
results. The first one follows from the identification of the 
most singular behavior in the temperature dependence of 
the various susceptibilities, as shown in Fig. for values 
of t' ± below and above the threshold for superconductiv- 
ity. 

An alternative way to determine the nature of the or- 
dered phase is to look at the wave vector dependence 
of the renormalized vertex functions close to the diver- 
gence. In Fig. |3 we have plotted the SDW vertex func- 
tion r,g(— k' ± , kj_, kj_, — k±) in the k±,k'j_ plane. In the 
SDW regime, only processes involving particle-hole pairs 
at wave vector (— — k± ~ Q±q are found to be singu- 
lar. On the other hand, near the superconducting tran- 
sition temperature, we obtain the separable form 

rs(— fc^, k'j_, k±, — k±) oc cos k'j_ cos k± . 



8 




FIG. 3: -r s (-k'_ L ,k' x ,k±,-k±) in the SDW regime (t' x = 
0.04t±) close to the divergence (a), in the SCd regime (t'± = 
0.048t±) at an intermediate stage of the flow (b) and close to 
the divergence (c). The vertical scale is arbitrary. 



of the variables + and (k± — k' ± ). Neglecting the 
weak (k± — k' ± ) dependence, the result is, schematically, 

g s (-k' ± ,k' ± ,k ± ,-k ± ) (32) 
1 2 

= £ M-ir 

n= — j^ + l 

x (cos nk' ± cos nfcj_ — sin nk' ± sin nfc^) . 

\n\ corresponds to the distance between the chains where 
the two interacting electrons are located. In the SDW 
regime, Fig. Eta), g s has the form of a 8 peak, so that 
all a n > will be equal. In the superconducting regime, 
Fig.Etb), the peak is slightly enlarged, so that a n > will 
be some decreasing function of \n\. We thus see that, due 
to the spin fluctuations, the effective interaction between 
particles forming a singlet pair contains attractive contri- 
butions at all chain distances. (Note that, according to 
our definitions Eqs. I|29|l and (|3U|) . an "attractive" inter- 
action corresponds to positive g s ,t) The most important 
one is the nearest-neighbor chain one, cosfcj^ (SCd), fol- 
lowed by sin2/c x (SCg), cos3/c_l etc., as it can also be 
seen from the calculation of the related pairing suscep- 
tibilities, cf. Fig. EJb). Note that, according to the re- 
lations (|3I[) . gs also gives an attractive contribution to 
the triplet channel. However, all the three components of 
a - spin-one boson - SDW fluctuation contribute to the 
superconducting coupling in the singlet channel, whereas 
only one contributes to the triplet channel. Antiferro- 
magnetic fluctuations thus favor singlet pairing as com- 
pared to triplet pairing. On the other hand, according to 
Eq. (|31|l . strong charge fluctuations should in principle 
be able to change this tendency. 



A plot of the singlet interaction amplitude T s gives es- 
sentially the same picture. From Eq. I|27|l . it is clear that 
the most divergent vertex part £ in such a situation is 
the one proportional to cos k± . 

Recall that the vertex functions T we calculate in our 
RG scheme are equal to the effective interactions g e $' 
of a Wilsonian low energy effective theory. Quite early 
during the flow, the vertex associated to the spin density, 
Tg(kj±), and hence gs(kj±_), develops a peak structure 
similar to the one close to the SDW transition, but less 
pronounced, cf. Fig.Etb). There are thus important spin 
fluctuations at temperature or energy scales above the 
transition temperature for superconductivity, as is con- 
firmed by the behavior of the associated susceptibility, 
Fig. EJb) . The same peak structure appears in the effec- 
tive interaction between electrons forming a singlet (or 
triplet) pair, g s (t)(kj±). We can decompose the peak of 
g s (—k' ± ,k' ± ,kj_,—k±) at intermediate cutoff A in terms 



III. INTERCHAIN INTERACTIONS AND 
TRIPLET SUPERCONDUCTIVITY: 
INCOMMENSURATE CASE 

In this section, we examine the role of interchain in- 
teractions in the phase diagram of the extended quasi- ID 
electron gas model. We will first consider the influence of 
interchain backward scattering g± and forward scatter- 
ing g^ separately, before studying their combined effect. 
The intrachain interactions are kept fixed to their val- 
ues used in the previous section. Concerning the nesting 
quality, the overall picture remains the same. For weak 
deviations from perfect nesting, we find density-wave in- 
stabilities, superconductivity for more important devia- 
tions, and - discarding possible Kohn-Luttinger effects 
at low temperatures beyond the numerical accuracy of 
our RG calculations - a metallic phase when the nesting 
is deteriorated even more. 
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FIG. 4: Low temperature phases versus gi , keeping g 2 — 
and §3 — 0. Circles indicate a SDW phase, squares a CDW 
phase, triangles SCd (cosfcx) and crosses SC/ (r cosfcx). In 
the region below the dotted line, spin fluctuations dominate 
over charge fluctuations in the normal phase. 

By switching on g/ gradually on the positive scale, one 
observes that at small values of .g/, the SDW phase re- 
mains essentially unaffected, whereas the transition tem- 
perature for d-wave superconductivity is considerably re- 



duced by the presence of a finite g^ (Fig. [SJ, and the 
region where superconductivity is stable shrinks (Fig.|4|. 
For higher g^ (<?/ of the order of ^ of the intrachain 
backscattering <?i), d-wave superconductivity turns out 
to be no longer stable and an /-wave triplet supercon- 
ducting phase appears. In this parameter range, (2fcp, 7r) 
charge fluctuations are strongly enhanced by interchain 
backward scattering (Fig. EJ). However, spin fluctuations 
remain important, and in a sizeable region of the phase 
diagram, triplet superconductivity is preceded in temper- 
ature by dominant spin fluctuations in the normal state 
(Figs. 0]and[|JJ). For the values of intrachain interactions 
used in this section, and <?/ = 0.105, one finds a max- 
imum T C °(SC/) « O.OOlto. - 0.2 K, that is of the same 
order of magnitude as that of the d-wave case. T C (SC/) 
increases with the amplitude of <?/, and the supercon- 
ducting phase widens (Figs. 0] and [5J. Once triplet su- 
perconductivity occurs, g^ starts to affect the density 
wave phase: The SDW state is suppressed and replaced 
by a CDW. The values of 3/ for which SC/ and CDW 
phases first appear depend on the values of intrachain in- 
teractions and increase with the value of the ratio (71/32- 



I 

The origin of the /-wave SC and CDW phases can be understood by considering the contribution of the g^s to 
the (bare) scattering amplitudes in the singlet and triplet particle-particle channels, as well as in the charge and spin 
channels: 

9c( k '± +<7-L> fc -L - Qx,k' x ,kx) = -^cosqll +2g£cos(q± + k' x - k x ) , (33) 
gs(k' x +q±,k±- q±,k' x ,kx) = 2g% cos(q ± + k' x - k ± ) , 



gj-(-k'x,k'x,kx,-kx) = 2(-.g 1 L - g£) cosk'x coskx + 2(g^ - g£) sin/c^ sinfc_L , (34) 
g^-(-k'x,k'x,kx,-kx) = 2(3/ - g%) cosk' ± cosk± + 2(-g 1 L - g£) sink' ± smk± . 

From these equations, it can be easily seen that the interchain repulsion g/ contributes positively to g^j at g^ = 7r, 
and therefore induces CDW correlations with a phase difference of n between neighboring chains. In the Cooper 
channel, g^ favors triplet /-wave and singlet d xy -wave pairing, whereas its contribution to singlet d x 2_ y 2-wave and 
triplet p y -wa.ve pairing is negative. As for g^ , it tends to suppress both singlet and triplet pairings on nearest- neighbor 
chains. 

In addition to this 'direct' contribution there is also an indirect effect due to the exchange of density fluctuations. 
Upon renormalization, (2kp, tt) CDW correlations are enhanced beyond the level expected from a mean-field treatment 
of gi~ These CDW fluctuations enhance triplet /-wave pairing but suppress singlet pairing, whereas the SDW 
fluctuations are well known to favor singlet pairing (see Eqs. |3ip). Eq. (|32|l shows that the latter reinforce d x 2_ y 2- 
but suppress d xy -wave pairing. Regardless of the values of the interaction constants gi, gi and g% , the CDW and 
triplet /-wave phases always appear almost simultaneously when 3/ increases. This suggests that CDW fluctuations 
(rather than the direct effect of g± in the Cooper channel [Eqs. J31J]) provide the dominant driving force leading to 
/-wave superconductivity. 



We now consider the effect of interchain forward scat- 
tering g^ alone, setting g/ = 0. From Eqs. (|34|l . it can 
easily be seen that g^ contributes negatively to both </j- 
and g^~ and then favors the suppression of the nearest- 
neighbor-chain Cooper pairing induced by spin fluctu- 
ations. Our results (Figs. [7| and [SJ show that this is 



indeed the case. We have seen, however, in section III Dl 
that SDW fluctuations can generate a smaller yet present 
attractive interaction between electrons on next-nearest- 
neighbor chains, which is not affected by g% ■ It fol- 
lows that when d-wave superconductivity is sufficiently 
weakened by g^, it is replaced by g- wave singlet pairing 
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FIG. 5: Transition temperatures for different values of g\ , 
when §2 — and g-j = 0: 

(a) gi = (phase sequence SDW— >SCd, continuous line), 
0.090 (SDW^SCd, dotted line) and 0.120 (CDW^SC/, 
dashed-dotted line); 

(b) gi = 0.090 (SDW^SCd, dotted line) and 0.105 
(SDW(->SCd)->SC/, dashed line). Note that in the latter 
case, the SCd phase is extremely narrow. 



(A r (k±) oc r sin2£;^). This is shown in Fig. [7J 

It is worth noting that the instability of the normal 
state with respect to superconductivity with high angu- 
lar momentum pairing can be seen as a result of the 
Kohn-Luttinger effect originally predicted for isotropic 
metals ^ At variance with more isotropic systems, how- 
ever, the RG results show that for a quasi-lD metal, the 
transition temperature of high angular momentum su- 
perconducting phases, like SCg, remains experimentally 
accessible. We find T™§* g ~ 0.002i_L - 0.4 K. Finally, 
the SDW phase remains nearly unaffected by g^r ■ This 
can be understood from Eqs. I|33|> . which indicate that 
its contribution averages out over the Fermi surface in 
the particle-hole channel of the RG equations. 

We now consider the combined effect of g^ and g^ . 
Fig. shows the results for g± = g^, where all aforemen- 
tioned phases appear. In the presence of both g/ and g£~ , 
d-w&ve pairing is suppressed even faster than by <?/ or 
g% alone. The appearance of /-wave pairing is retarded 
by <72~ which, as previously mentioned, is detrimental 
to nearest neighbor chain pairing (Eq. ijslljl). However, 
for our choice of intrachain interactions, no triplet phase 
with pairing on next-nearest-neighbor chains is found for 
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FIG. 6: Temperature dependence of the susceptibilities in 
the normal phase above the SCd phase ((a) g x — and (b) 
<7i = 0.08) and above the SC/ phase ((c) gi — 0.10), keeping 
92 — 0, <73 = and t' ± = 0.056tx- The continuous line 
corresponds to SDW, the dotted line to CDW, the dashed 
line to SCd and the dashed-dotted line to SC/ correlations. 



<7i = i since g^~ favors nearest-neighbor triplet pairing 
at the outset. As for CDW's, they are found to occur 
at slightly higher values of g^~ in the presence of a finite 
02 • As mentioned before, g^ has no important effect in 
the particle-hole channel alone. The slight suppression 
of the CDW due to g^; must therefore come from the ID 
regime, where the correlation channels are coupled. This 
can be checked from the RG equations for t± = given 
in appendix 1X1 

The phase diagram (Fig.[5J) depends quantitatively and 
qualitatively on the bare intrachain interactions. When 
the ratio g\j gi increases, higher values of the interchain 
interactions are necessary for obtaining CDW and SC/ 
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FIG. 7: Low temperature phases for g 2 varying, keeping g± = 
and §3 = 0. Circles indicate a SDW phase, triangles SCd 
and diamonds SCg (r sin2fcx). 
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FIG. 8: Transition temperatures for different values of g 2 , if 
§i — and §3 = 0: §2=0 (SDW— >SCd, continuous line), 
0.16 (SDW^SCd, dashed line) and 0.20 (SDW^SCg, dotted 
line). 



phases because Tc is more negative at the outset. For 
.<?i = .92, we even find a triplet SCf order parameter of the 
form A r (fcx) sm2k± 1 corresponding to second- nearest- 
neighbor chain pairing, instead of A r (k±) oc rcosk±. 
This is coherent with the fact that, in the ID regime, the 
renormalization of T^ is proportional to r^Tc, so that a 
more negative Tc reduces and hence nearest-neighbor 
chain triplet pairing (appendix ^ . 



IV. EFFECT OF UMKLAPP PROCESSES 

Conductors like the Bechgaard and Fabre salts are 
slightly dimerized in the direction of the organic chains. 
It follows that at low energy or temperature, the hole 
band can be considered as effectively half filled rather 
than quarter-filled and this gives rise to Umklapp scat- 
tering processes with amplitudes gs and g$. To leading 
order, the bare amplitude of 173 ~ giAjj/Ep is propor- 
tional to the dimerization gap Ad ^LSL which yields a 
g3 that is rather weak as compared to gi.2- We assume 
a similar ratio between the interchain couplings g$ and 



FIG. 9: Low temperature phases for gi = g 2 , (?3 = 0. Circles 
indicate a SDW phase, black squares a CDW phase, triangles 
SCd (cosfc^), white diamonds SCg (r sin2fcx) and crosses 
SC/ (r cos Below the dotted line, spin fluctuations dom- 
inate over charge fluctuations in the normal phase. 
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Intrachain Umklapp processes enhance the formation 
of 'site' centered SDW correlations and 'bond' centered 
CDW correlations and weaken bond SDW and site CDW 
fluctuations, as can be seen from Eqs. Q25[l. Thus, con- 
trary to the incommensurate case, site and bond density- 
wave orders no longer join together and must be consid- 
ered separately with their distinct strengths at half filling. 
In addition, intrachain Umklapp scattering reinforces the 
effective interactions gc and gs at momentum transfers 
(2k p, tt), thus increasing spin and charge density fluctua- 
tions of both site and bond type. These effects are known 
in the ID case (see e. g. Ref . IsqI and also Appendix 0, 
as well as from calculations restricted to the particle-hole 
channels (Ref. Ffol and Appendix 0. Even though Umk- 
lapp processes do not enter the particle-particle channel 
directly, they do enhance Cooper pairing, because they 
amplify the peak at (2k f, t) in the effective interactions 
gs.c an d thus the attraction between electrons on neigh- 
boring chains [see Eq. (|32|l ]. Since density- wave correla- 
tions are enhanced with respect to the incommensurate 
case, transition temperatures are also higher^ This is 
true for the formation of density-wave states as well as 
for superconductivity (Fig. \M&)). The ratio T c °/T c m s °^ 
remains essentially unaffected by g%. The critical t' ± 
needed to destroy the density-wave phase is also higher 
in the presence of Umklapp scattering. 

Let us now turn to the effect of interchain interactions 
in the presence of Umklapp processes. We will first con- 
sider the effect of a finite intrachain 173, for g^ = > 
and g% — 0. A comparison of Figs. 151 and lllf a) shows 
that the regions of d-wave and /-wave superconductivity 
are now enlarged, and there is no more SCg phase corre- 
sponding to next-nearest-neighbor chain singlet pairing. 
This is a consequence of the fact that density- wave corre- 
lations are reinforced by Umklapp scattering. A stronger 
#2 would then be needed to destroy nearest-neighbor 
chain pairing, but owing to the condition g^ = g^ the 
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FIG. 10: Transition temperatures when Umklapp processes 
are taken into account. 

(a) Without interchain Umklapp scattering: §3 = 0, gf = 
(SDW— >SCd, continuous line), 

§3 = 0.02, gf = (SDW^SCd, dashed line), 

§3 = 0.02, gf = gf = 0.1, gf = 0, (SDW->SC/, dotted line). 

(b) Effect of interchain Umklapp processes: §3 = 0.02, gf — 
gf = 0.1, <?/ = (SDW^SC/, dotted line) and 

gf = ^~gf (SDW->SCd->SC/, continuous line). 



SC/ phase is also reinforced and in turn stabilized once 
the SCd phase is suppressed. 

When we finally add interchain Umklapp scattering 
gf, the picture does not change significantly (Fig. lTl7 b)). 
The only effect is that the occurrence of SC/ and CDW 
phases takes place at slightly higher values of the inter- 
chain interactions in the phase diagram, whereas transi- 
tion temperatures are scarcely lower (Fief. HUf b)). In the 
ID regime (A » t±), gf enhances gc (appendix [XJ, 
but this is only a second order effect. On the other 
hand, the bare contribution of gf, Eq. ||BJ|, reduces the 
peak in the effective interaction ga(k[±, k' 2±1 /c2_l, fcix) al 
k[j_ — ki± = it, which is connected to the charge den- 
sity fluctuations (appendix E) . This influence is linear 
in gf and therefore stronger than the ID effect, thus ex- 
plaining the overall weakening of CDW and, in turn, of 
SC/ correlations. There is no direct effect of gf on SDW 
correlations. However, the RG equations (Appendix 
show that gf weakens the renormalization of #3 in the 
ID regime and in turn the increase of transition temper- 
atures in comparison to the incommensurate situation. 

In Fig. El we show the typical behavior of the most 



FIG. 11: Low temperature phases in the presence of Umk- 
lapp processes, with gf = (a) and gf/gs = gf/gi (b), 
taking gf = gf. Circles indicate a site-SDW phase, squares 
a bond-CDW phase, triangles SCd (cosfcx) and crosses SC/ 
(r cos k±). Below the dotted lines, spin fluctuations dominate 
over charge fluctuations in the normal phase. 



important susceptibilities as a function of temperature 
in the normal phase above the superconducting phases. 
Fisf. ll2f a) shows that for sufficiently important interchain 
interactions, triplet correlations are already strongly en- 
hanced, although the ordered phase still corresponds to 
spin singlet pairing. SDW correlations are always impor- 
tant. They are the dominant fluctuations in the normal 
phase in the major part of the parameter range we have 
explored, sometimes even above the triplet SC phase (cf. 
Figs.E2Ic) and ITI) . 



V. DISCUSSION AND CONCLUSION 

In this work, we have determined the possible elec- 
tronic phases of the extended quasi- ID electron gas 
model that includes both intrachain and interchain re- 
pulsive interactions, interchain electron hopping and the 
influence of nesting deviations. Our results reveal that 
in correlated quasi- ID metals both for zero and non zero 
Umklapp scattering, interchain interactions can act as 
a key factor in expanding the range of possibilities of 
ordered states compared to the case where only repul- 
sive intrachain interactions are present ^2*22i2i At large 
momentum transfer, the interchain electron-electron cou- 
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FIG. 12: Temperature dependence of the susceptibilities in 
the normal phase above the SCd phase ((a) t'± = 0.152tx, 
gi = 0.08) and above the SCf phase ((b) t' ± = 0.152t x and 
(c) t' ± = 0.1761l, both for gi = 0.12), taking g 3 = 0.02 



and g 2 = g-y 



l -g s . The continuous line corresponds to 



SDW, the dotted line to CDW, the dashed line to SCd and 
the dashed-dotted line to SCf correlations. 



pling acts as a short-range interaction that is responsible 
for the enhancement of CDW correlations, consistently 
with what was found long ago in the absence of inter- 
chain hopping4iaa For finite fj_, however, interchain in- 
teraction leaves the amplitude of SDW correlations es- 
sentially unaffected and a relatively small critical value 
of gi coupling is then needed to make the CDW ordered 
state possible, and this, on equal footing with the SDW 
phase, which is known to dominate the phase diagram 
when only repulsive intrachain interactions are present. 
When nesting deviations are cranked up beyond some 
threshold value t\ , density- wave order is suppressed and 



interchain interactions turn out to affect correlations of 
the Cooper channel too. Thus an important conclu- 
sion that emerges from this work is the gradual suppres- 
sion of interchain d-wave pairing when repulsive short- 
range interchain interaction increases. As a result of the 
growth of CDW correlations in the normal state, it ul- 
timately yields the stabilization of a triplet SC/ super- 
conducting phase corresponding to an order parameter 
A r (fcj_) = rAcosfcj_ having nodes on the warped Fermi 
surface that are at the same locus as for the SCd case. 
The normal phase is still dominated by strong SDW fluc- 
tuations over a sizable region of the phase diagram in this 
sector. 

It is interesting to consider how far the RG results of 
this work are applicable to quasi- ID organic conductors. 
In the case of the Bechgaard (TMTSF) 2 X salts with cen- 
trosymmetric anions X for example, the observation of 
a SDW-CDW coexistence below the critical pressure for 
superconductivity22i£2i£i indicates that interchain elec- 
tron repulsion is a relevant interaction in these materials 
besides intrachain interactions and interchain hopping. 
Accordingly, the RG calculations show that a relatively 
small and realistic amplitude of repulsive interchain in- 
teraction is sufficient to bring the stability of CDW or- 
der close to SDW, indicating that this part of interaction 
would indeed play an important role in the emergence 
of density- wave order in these materials. The suppres- 
sion of the SDW state (and presumably of CDW as well) 
followed by the emergence of superconductivity is well 
known to constitute the closing sequence of transitions 
that characterizes virtually all members of the Bechgaard 
and Fabre salts series as one moves along the pressure 
scale. Since pressure introduces alterations of the elec- 
tron spectrum, it prompts deviations from perfect nest- 
ing. In our model, these simulate the main effect of pres- 
sure, which together with a reasonable set of parameters, 
yield a 'pressure' profile of the critical temperature that 
agrees quite well with the characteristic variation seen in 
experiments^ As to the nature of the superconductivity 
in these materials, our results based on a purely elec- 
tronic model, indicate that given the observation of the 
close proximity between SDW and CDW ordered states, 
not only SCd but also triplet SC/ order parameter be- 
come serious candidates for the description of the su- 
perconducting phase in these compounds (Figs. 14191 and 
lll|l . However, as pressure also affects the normalized 
amplitudes through the band width and the dimer- 
ization gap (Umklapp), and gr- through the interchain 
distanceji^ the actual trajectory in the phase diagram 
under pressure cannot be determined with great preci- 
sion. It follows that besides the possibilities SDW— > SCd 
or SDW — > SC/, sequence of transitions such as SDW — > 
SCd — > SC/, where one can pass from singlet to triplet 
SC order under pressure, cannot be excluded. It is worth 
remarking that in this sector of the phase diagram, the 
addition of a small magnetic field - as actually used in 
many experiments^^ - or accounting for the small but 
yet finite spin anisotropy would tend to tip the balance 
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in favor of a triplet order parameter £L£L 

Experimental features of the normal phase also argue 
in favor of this region of the phase diagram for the Bech- 
gaard salts. This is the case of the puzzling growth of 
CDW correlations seen in optical conductivity in the low 
temperature part of the metallic phase above the super- 
conducting transition^ CDW correlations are found to 
be significantly enhanced in a temperature region where 
NMR experiments reveal the existence of strong SDW 
correlationspiii^ This feature cannot be captured for re- 
alistic intrachain interactions alone. It can find, however, 
a natural explanation in the framework of the extended 
quasi- ID electron gas model, for which interchain inter- 
actions can boost the amplitude of CDW correlations be- 
sides those of the SDW channel that are kept essentially 
unchanged. 

In the case of non-centrosymmetric anions (e. g. 
X=C104), our approach should be refined in order to take 
into account the doubling of the unit cell in the transverse 
direction due to the anion ordering taking place below 
24 K in the normal phase. The concomittant reduction 
of the Brillouin zone yields two electronic bands at the 
Fermi level and in turn multiple nesting vectors^2& An 
accurate description of this Fermi surface should be in- 
corporated in the RG approach, and a modification of 
the node structure for the SC gap is expected^ 7 - A previ- 
ous (simplified) RPA-like calculation has shown that the 
nodes of a d-wave SC order parameter, which are located 
at k± = ±7r/2, are precisely found where a gap opens 
due to CIO4 anion ordering, thus making the SC phase 
effectively nodeless at low temperature S 

The results presented in this work may also be relevant 
for the phase diagram of other series of quasi- ID organic 
superconductors. In this matter, the case of the two- 
chain compounds TTF [M(dmit) 2 ] 2 ( M=Ni, Pd) is of 
particular interest. These compounds are characterized 
by an incommensurate CDW state that takes place on 
the M(dmit)2 stacks at low pressured At some critical 
pressure, they become superconducting. 79 *^, The varia- 
tion of the critical temperature under pressure has been 
analyzed in detail in the case of M=Pd, which shows 
close similarity with the one found for the Bechgaard 
and Fabre salts The temperature scale for the CDW 
instability at low pressure is relatively large, however, 
and owing to the pronounced anisotropy of the band pa- 
rameters in these compounds* 7 ^! this indicates that the 
interchain interaction is likely to be a key coupling in 
the stabilization of a 3D ordered state in these materi- 
als. According to our model, these conditions would be 
favorable to the existence of a triplet SC/ state in these 
systems under pressure (see e.g. Fig.[5J). 

Another result that is highlighted by our analysis in the 
incommensurate case is the occurrence of a singlet SCg 
state when the interchain electron scattering dominates 
for small momentum transfer. The 5^ coupling tends 
to suppress electron pairing between the first nearest- 
neighbor chains, which therefore suppresses SCd or SC/ 
type of superconductivity. However, an instability of the 
normal state remains possible. It results from the oscil- 



lating tail of density-wave correlations in the transverse 
direction which favors longer range pairing between elec- 
trons separated by more than one interchain distance. 
The fact that an instability of the normal state persists 
in the Cooper channel illustrates how the quasi- ID ge- 
ometry for electrons, with its inherent interference be- 
tween Peierls and Cooper channels, is prone to magnify 
the Kohn-Luttinger mechanism for the (Cooper) instabil- 
ity of a Fermi liquid for repulsive interactions^ 7 - One can 
easily infer from our results that the addition of longer 
range interchain interactions will frustrate short-range in- 
terchain pairing and shift it to larger interchain distances, 
thus unfolding possibilities of superconductivity at even 
larger angular momentum pairing. 
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APPENDIX A: RG EQUATIONS IN LIMITING 
CASES 

In order to obtain the one-dimensional limit of the RG 
equations given in section lll Bl we neglect interchain hop- 
ping (t_L = t' ± = 0) and take into account intrachain 
and nearest-neighbor-chain interactions only. We thus 
obtain^ 

= _ ( f(^)) 2 _2[(f^) 2 + (f^) 2 ] , (Al) 

rr = 4(fr) 2 +i(fn 2 , 

r<">> = f( 1D) [-f( 1Z5) + 2f( 1Z5) ]-4f^, 
H = f$[-2rw+fw + fi]-f?°)fi. 

Alternatively, we may restrict the renormalization to 
a given correlation channel in certain situations. For 
the particle-hole channel, it is convenient to rewrite the 
RG equations in terms of the vertex functions Tc,s an d 
T^' S \ as defined in eqs. (|29|l . Remember however that 
Tg and are not independent: 
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(a) 




FIG. 13: Typical pictures of the interactions 
T s (-k' ± ,k' ± ,k±,-k±) (a) and T 3 {-k' ± ,k' ± ,k±,-k±) (b), for 
initial values gs = gc = 0.2, gz — 0.02, if the renormalization 
is restricted to the Peierls channel. The result for Tc is 
identical to that for Tg. 



irVix, k' 2 ±,k 2 ±,k 1± ) = -2T\? } (k' 1± ,k 2± ,k 1± , k 21 _) + T { 3 b) (k' 1± ,k' 2± ,k 2± , k 1± ) . 



We obtain 

fc{k' 1 j_ 7 k' 2±7 k 2 ±,k 1 j_)\ p = -^-^2B v (k ± ,k' 1± - k 2 ±) 

x {f c (k' 1A _,k ± ,k 2± , k' ± )t c {k' ± ,k 2± ,k±,k 1 j_) 
+fP (k' 1± ,k ± ,k 2± , k' ± )fP (k' ± , k' 2± , k± , k 1± ) } , 



^s(k' 1 _ L ,k' 2± ,k 2 ±,ki±)\ v = — ^B v {k^,k' 1 ^-k 21 _) 

^ k± 

x {fs(fcij_, k±, k 2 j_, k'j_)t s (k'j_, k' 2± , k±, ki±) 
+Ti s) (k' 1± , k ± , k 2± , k' ± )Ti s) (k' ± , k' 2± , k ± , k 1± ) } , 



1G 



?l s \k[ x ,k' 2± ,k 2X ,k lx )\ v = J-J2 B r( k ^ k 'i±- k ^) ( A4 ) 



x 



N± 

K± 

l - x {-t c (k' 1± ,k ± ,k lx ,k' ± )T^ ) (k' ± X±,k±,k 2 ±) 

-^3 C) (K±> fc -L> fc ix, k' ± )t c(k'j_, k' 2± , fc_L, fc 2 x) 
+f s (fc' lx , fcx, fci ±1 *i )f I s } (fci , fc 2± , k x , k 2± ) 
+ff ) {k[ ± ,k ± ,k lx ,k' ± )f s (k' ± , k' 2± ,k x ,k 21 _)} 

x {f s (fei ±I fe 2 x, fc ± )f < s } (fc x , fc 2X , k x , feix) 
+f ^ (fci x , fcx, fc a ±, fc x )f s (k' x , k' 2X , k ± , fcix) } • 

These equations are remarquably symmetric with respect to spin and charge density correlations, cf. e.g. the results 
when we take gs — gc from the beginning (Fig. I13jl . For low temperature and energy cutoff A, the particle-hole loop 
integral B-p(q±) is strongly peaked for momentum transfers equal to the best nesting vector, i. e. q± = n. As can 
be seen from the preceding equations, this generates a peak at k' 1± — k 2 x — n in T$ (responsible for the SDW) and 
in Tc (responsible for the CDW). In T^, peaks are created at k' 1± — k 2 x = tt as well as k' 1± — k±x = 7r. The peak in 
T% for k[ x — k 2 x = 7r enhances the spin density correlations, whereas the one for k' ix — fcix = 7r supports the charge 
density correlations. 

Let us now consider the influence of interchain interactions [Eq. |JBJ|] on these structures. Interchain Umklapp 
scattering g^- reduces the latter peak and thus weakens the charge fluctuations, whereas it has no direct influence on 
the spin correlations, g^ and g 2 have no direct effect on the spin fluctuations either. Interchain backward scattering 
<7i , on the contrary, reinforces the peak in gc at k' 1± — k 2 ± — tt and thus supports the formation of charge density 
waves. 

We finally give the RG equations restricted to the particle-particle channel. They are most conveniently written in 
terms of the singlet and triplet pair interactions, a = s, t: 

r a (fc ±1 /c ±2 /cx2fcxi)| c = ^^5 c (fcx,<zcx) ( A5 ) 

fc_L 

xf a (k' ±1 k' ±2 kxk' ± )T a (k±k' x kx 1 kx2) . 

k ± = -k ± +q c± 
1C±= k ±l+ k ±2 
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